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ABSTRACT ; 

Procedures for the practical implementation of the hybrid coordinate 
methods of dynamic analysis of flexible spacecraft in application to 
vehicles of realistic complexity are briefly documented, with supporting 
examples. 
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INTRODUCTION ; 

Reference 1 provides a detailed derivation of equations representing 
the dynamic behavior of spacecraft idealized as collections of rigid bodies 
and flexible appendages; these equations employ a hybrid system of co- 
ordinates, including discrete coordinates for rigid bodies and distributed 
or modal coordinates for elastic appendages. The objective of Reference 1 
is the development and dissemination of a body of theory, with only minimal 
consideration being given to the ways in which this material might be used. 
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In the period following the preparation of Reference 1, serious 
attention has been given at JPL to the development of the practical 
utility of the hybrid coordinate method by exploration of the various 
ways in which it can be used in the practical design of attitude control 
systems for flexible spacecraft. The application of hybrid coordinate 
design procedures to specific spacecraft has been described m separate 
papers as these methods have evolved (References 2,3); the present report is a 
compilation and documentation of the present status of this continuing 
process of developing and refining design tools which have become available 
only with the introduction of the hybrid coordinate method of Reference 1. 

It may be remarked that the hybrid coordinate method has been widely 
adopted outside of JPL, and in project applications design procediares have 
evolved in other organizations. This Memorandum is however largely con- 
fined to those design methods developed and applied within JPL. 

Detailed technical material is reserved in this Memorandum for three 
appendices, each treating a specific example. The emphasis in the body 
of this report is on uhe qualitative description of design procedures 
and the comprehensive strategy of approaching the task of attitude control 
system design for spacecraft with large flexible components. 

DESIGN PROCEDURES ; 

A feature of the hybrid coordinate formulation of equations of motion 
which distinguishes them from those obtained by the vehicle normal mode 
method is the retention in the former of explicit attitude variables for 
that rigid portion of the vehicle for which attitude control must be 
maintained. Because these variables (perhaps a set of three attitude 
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angles) remain explicitly in the equations, it is a simple matter- "when 

the equations are linearized to obtain a transfer functxon relating the 

transforms of these attitude variables to transformed control system torques, 

autopilot gimbal angles, or whatever measure of control system actuation 

is most convenient. The influence of vehicle flexibility, as reflected 

in the response of any number of flexible appendage modal coordinates, is 

manifested in the transfer functions for the attitude variables . Thus the 

hybrid coordinate method is particularly well-adapted to the many control 

* 

system design procedures which employ transfer functions. 

At JPL, transfer functions have in the past been used primarily in 
the form of scalar (single-axis) transfer functions for the generation of 
root locus plots, A computer program developed by E. Kopf and R. Mankovitz 
permits rapid and inexpensive construction of root locus plots, which have 
proven to be extremely useful for preliminary design of flexible vehicle 
control systems (see Appendices A, B, C). 

'Other organizations (e.g. Hughes Aircraft Company) have used Reference 1 
to obtain transfer functions which they have employed in generating phase- 
gain plots, and any of the other conventional alternatives to root-locus 
plots can equally well be utilized. 

For a rigid body, the assumption of uncoupled single-axis response is 


This is not to say that transfer functions cannot be obtained from equations 
employing either exclusively discrete coordinates or exclusively modal co- 
ordinates, since any linear set of equations will suffice, if one requires 
transfer functions only for the vehicle attitude variables however, with 
deformation variables occupying a secondary importance, then transfer 
functions obtained from hybrid coordinate equations are both conceptually 
simpler and easier to calculate. 
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jeopardized only by the existence of products of inertia for the control 
axes or the presence of stored angular momentum normal to the control axis. 
For a flexible s/C, one must add to these considerations the likelihood 
that any given mode of vibration of a flexible appendage -will also cause 
cross -axis coupling. Except for the simplest of idealized S/c (such as 
that in Appendix A), it must be expected that single -axis response studies 
for flexible s/C will be less valid in predicting vehicle response than 
has been the case with the nominally rigid S/c of the past. Whereas it 
has been a common past practice in S/c control system design to proceed 
directly from linearized, single-axis response studies using root locus 
plots to digital simulation of coupled, nonlinear equations, it is recom- 
mended that an interraediate step be inserted when vehicle flexibility is 
expected to introduce cross-axis coupling. This step involves the analysis 
of coupled equations for which the approximations of linearization have 
been retained. 

The use of transfer functions does not preclude the incorporation of 
cross-axis coupling, but such coupling does complicate their use as a 
design tool. One can select one of the three axes (say axis 1 ), and 
establish a control law for the torque T^^ about this axis to maintain 
control of the corresponding attitude angle and then proceed to the 
second axis and devise a control law for Tg with the recognition that 
©2 is influenced not only by Tg but also by T^, which is introduced ac- 
cording to the control law as a consequence of the ©^ induced by Tg. By 
the time one faces the prospect of choosing a control law for T^ to control 
©2, the influences on ©^ of (available as a function of ©^) and Tg 
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(available as a function of 0^) greatly complicate the design task, and 
the designer may wonder if his job would have been easier had he started 
with the 3“3^is; but there are no conceptual obstacles to this approach. 

At JPL the problem of analyzing coupled linear equations has been 
confronted directly with a digital computer program for the determination 
of eigenvalues and eigenvectors, eschewing coupled-axis transfer functions. 
After the investment of some effort in obtaining a suitable eigenvalue- 
eigenvector program and acquiring the necessary skills in its use , this 
has become a useful design tool for flexible vehicle attitude control 
systems. In Appendix C the results of eigenvalue studies for the TOPS 
vehicle are presented, and compared with preliminary conclusions based on 
single-axis response root locus studies. 

The final analytical and computational step in control system design 
is a computer simulation of the systan of coupled, nonlinear equations of 
motion; at JPL this stqp has been taken via digital computer numerical 
integration (see Appendices B and C). 

Thus at JPL the design of attitude control systems for complex flexible 
spacecraft has evolved into a three-stage process: (1) Linearized single- 

axis response studies using computer -generated root locus plots, (2) Linear- 
ized coupled system equation studies using an eigenvalue -eigenvector 
computer program, and (3) Nonlinear coupled system simulations via digital 
computer. 

The program EVPEVT with error bound print-out was finally adopted, with 
the pracrice of normalizing variables to improve the conditioning of the 
input matrix until acceptable error-bounds were attained. 


JPL Technical Memorandum 33-493 



In each stage of analysis, appendage flexibility is represented in 

the hybrid coordinate approach by means of modal coordinates, and this 

woiild appear to require an extensive modal analysis of a.n flexible 

structure. This is a substantial mathematical modeling and computing 

task, involving the personnel and computer programs of a structures group. 

In the early phases of control system design, the structural design is 

generally also in a preliminary phase, and detailed modeling is um^arranted. 

This obstacle to efficient system design has been neatly circumvented by 

the development of an appreciation of the physical significance of certain 

matrices which embody the influence of appendage flexibility on control 

system design. Specifically, it has been noted that the influence of a 

flexible appendage in a given mode of vibration on control system operation 

can be characterized by a frequency o, a damping ratio and the contribution 

T 

of the mode to a matrix defined in Reference 1 at 6 6 . This matrix has 

the units of moment of inertia, and represents (in a certain sense) the 

contribution of the flexible appendage to the vehicle inertia matrix. 

Although one must have a rough appreciation of the mode shape for a 

particular mode of vibration in order to estimate the contribution of that 
T 

mode to 5 6 , this estimate can be made in many cases without actually 

performing the nodal analysis. For preliminary design pruposes it is there- 

T 

fore often possible simply to guess values of a, ^ 6 for an 

hypothesized mode of appendage vibration, >/hile such a guess will often 
require substantial revision as the vehicle design firms up and digital 
computer modal analysis becomes appropriate, this alternative is apt to 
be far preferable to the old practice of basing preliminary design on the 
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idealization of the vehicle as a rigid body. 

As the design process converges, and fairly accurate characteristics 
of control system and structure can be ascertained, it becomes increasingly 
important _ that simulations incorporate those modes of appendage vibration 
which influence the control system. At the same time, there is ari 5 )le 
motivation to exclude as many modes of vibration as possible from the 
simulation, due- in general to the cost in computer running time of each 
additional coordinate, and in particular to the inordinately high cost in 
running time and computer accuracy of retaining high frequency nodes of 
vibration. 

In addition to the usual procedure of truncating modal coordinates 
by means of a frequency criterion, retaining only those modes with frequencies 
in the control system band -width, two more techniques have been developed 
at JPL for guiding truncation. 

T 

As noted previously, the contribution of a given mode to 6 8 is a 
measure of its influence on the attitude control system. This matrix 
should therefore be computed and printed out as part of the modal analysis 
task undertaken by the structural mechanics experts. It then becomes 
possible simply by inspection of the elements of a 3 x 3 matrix for each 
mode to determine the extent of its influence on control about a given 
control axis, and the degree to which it causes cross-axis coupling. 

The final criterion for modal coordinate truncation in advance of 
digital computer numerical integration is the influence of an individual 
mode on system eigenvalues. It is much cheaper to explore the consequences 
of various altemative coordinate truncations with an eigenvalue program 
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than with a numerical integration routine. 

It should be noted explicitly in this Memorandum that the hybrid 

coordinate method relies heavily upon the cooperative interaction of 

personnel in control systems and structural mechanics organizations. Wot 

only are the structures people and their computer programs essential to 

the detailed task of appendage modal analysis, they also are invaluable 

T 

in helping to foimulate the rough estimates of ct, and 5 6 that go into 
preliminary control system design. At JPL there is a strong emphasis on 
promoting technical communication at the working level between individuals 
in Structural Mechanics (Division 35) and Guidance and Control (Division 3^+) 
and it is very important that this kind of joint effort be encouraged else- 
where. Larger aerospace organizations have on occasion been handicapped 
in their attempts to exploit the advantages of l^rbrid coordinate analysis 
by apparent difficulties in fostering continuing cooperative technical 
interaction between structures and controls groups, 

APPEIIDIX A — AW IDEALIZED TEST VEHICLS : 

In order that the development of methods of control system design not 
be obscured by the mechanical complexity of the object being controlled, 
a simple, idealized test vehicle (Figure 1) has been devised to permit 
physically meaningful interpretation of all phases of the analysis. Figure 1 
sho■^^s four elastic members attached to a rigid body B which may also con- 
tain a rigid, symmetric rotor. The elastic bodies are massless except 
for a tip mass on each, normalized to unity (weight .454Kg) The central 
body B (including its rotor) weighs (L3Kg); so the vehicle mass is distributed 
equally between the rigid part and the four flexible attachments, which 
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Figure 1. Idealized test vehicle 

together comprise "what will be called a single appendage. The tip masses 
are particles, and the central body is inertially spherical, with weight 
moment of inertia (.758 and dimensions small relative to beam 

lengths L and A, chosen as(l.2m) and (.6m), respectively. Beam stiffnesses 
are chosen to make the long beams have natural frequencies of O.90 and I.81 Hz, 
respectively, in the 1-2 and 1-3 planes of the vehicle when attached to 
a stationary base B. The corresponding "cantilever mode" natural frequencies 
of the short beams are 2.71 and 4.52 Hz in the 1-2 and 2-3 planes, re- 
spectively. For the investigations reported in this appendix, all beams 
are assumed longitudinally inextens ible . 

The raodal coordinates in the matrix if of Eqs, (224) and (216) of 
Reference 1 are not the cantilever modal coordinates of the individual 
beams, nor are they normal -mode coordinates of the total vehicle. (The 
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presence of the rotor ntakes the vehicle normal-modes complex, and these 
modes are not used here in order to confine the analysis to real numbers.) 
The coordinates, in 7] establish the response of the appendage in modes in 

V 

which they A'fould vibrate freely if the rigid body was free in translation 
but constrained against rotation. V7ith this interpretation, it becomes 
evident that there must be eight coordinates in the ’column matrix T), and 
for this simple system one can almosi. guess the mode shapes. 

Using established digital computer eigenvalue -eigenvector routines, 

one can determine the eight natural frequencies and mode shapes quite 
* 

accurately. Results are presented schematically in Figure 2, and queintita- 
tively in Table 1, In the latter, node manbers and numbers associated 
ifith "d.o.f." (degree of freedom) are to be obtained from the first: diagram 
of Figure 2, 

For purposes of preliminary design, it may be desirable to make 
exploratory transfer function calculations before the appendage is fully 
defined and subjected to modal analysis. Examination of the matrix transfer 
function in Eq. (292) of Reference 1 reveals that of the appendage 
properties only the matrices a, and 5 are required. The autopilot 
transfer function does involve other appendage properties (see Reference 3 ), 
but in the first approximation Equation ( 292 ) VJill suffice for autopilot 
control also. Since ^ is actually assigned and not calculated, only a 
and 5 must be calculated for application of the transfer function in 
this equation. 

Appreciation is expressed to Mr, John Garba of JPL, Division 35 , who 
actually performed these computer operations, using the SAMIS program. 
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For this simple system one can readily estimate ^ 0.9 Hz, since 
this is the input cantilever mode natural frequency of the longest beams , 

-JT_ 

Furthermore, one can utilize the- physical interpretation of 6 5 offered 

in the previous section, and guess that in the first mode only the long 

beams will participate (as in Figure 2); therefore be given 

2 ^ 

by 2(.U54Kg)(l.2m) ^ (l.35Kg -m2), since this figure represents the con- 

tribution to of the tv;o long beams. One might for this simple system 
easily guess as well the third, fifth, and seventh mode frequencies and 
contributions to 6 since these are asypnetric modes which have 
frequencies corresponding to the known cantilever mode frequencies. The 
physical interpretation of 6 6 tells us that the symmetric modes [modes 
2, 4, 6, and 8 in Figure 2) contribute nothing to this matrix, so they can 
be completely ignored. 

T 

Confirmation of these estimates for frequencies and 6^ 6^ contributions 
comes from the computer-generated value of 6j using the definition in 
Eq. ( 287 ) of Reference 1. The results for 5 and 6 6 (with no truncation) 
are given by 

0 0 -2a b 0 0 

0,0 0 6^6 = 0 4 b 0 

0 2a 0 0 0 5b 

0 0 0 , 

6=0 0 a 

0 0 0 

a 0 0 

0 0 0 I 
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■where a = [.01 Kg-n"^); b = (.3^- Kg~Jn^); and computer generated numbers 

below (353f>9)have been replaced by zeros. (For unit consistency it must be 
T — T -1 

noted that 6 6 is actually 5 E 6, with E a unit inertia matrix.) For this 

simple system, the estimated values are entirely correct. (The predicted 

^ 3^3 checks (-I. 16 ) squared, and the calculated 5^5 is 

precisely the contribution of the appendage masses to the vehicle inertia 

matrix. Furthermore, the even-numbered rows of 6 are null, so these m.odes 

T 

do not contribute to 6 6 .) 

Having obtained, by estimation or calculation, the necessary param- 
eters of the flexible vehicle, the control system designer can begin the 
sequence of response analyses required for selection of a control system. 

In general, the first step might be the rapid construction of a number of 
root-locus plots for single-axis response, the second step might then be 
a series of eigenvalue studies for a smaller range of prelimina.ry control 
system designs, and the final step would be a numerical integration of the 
most general system of equations of motion, to confirm the final control 
system design. For the simple test vehicle treated in this section, how- 
ever, each appendage modal vibration contributes to the vehicle response 
about one axis only, so the three scalar equations of vehicle motion 
(Eq, 224 of Reference 1) are coupled only by the "gyroscopic coupling" of 

A 

the rotor. Furthermore, response sbout the b^ axis (See Figure 1) 

A A 

is completely uncoupled from responses 0^ and ©^ about axes b^^ and 
so root locus plots for this axis must indicate precisely the same response 
that would be obtained from eigenvalue calculations or numerical integrations. 
The puzposes of this paper are best served by restricting the test vehicle 
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analysis results to the presentation of root locus plots for^G^* In the 
next two appendices eigenvalue analyses -and integration results are 
presented for space vehicles of realistic complexity. 

Figures 3-6 portray root-locus plots for G^ response of the test vehicle, 
under various assumptions. The dashed-line loci on Figures 3-5 sire based 
on the assumption that the vehicle is rigid , so they provide for three 
simple control systems a basis for evaluating the influence of flexibility. 
Figures 3-6 introduce flexibility in the first mode of vibration (see 
Figure 2 and Table 1), In this mode the normalized reduced inertia 
E 2 is 1 - (6263/I3) = 1 - [l.35/( .76 +1.68)] = 0.45. The natural 
frequency a is O.90 Hz, although in the plots this value is normalized 
to unity. The damping ratio is assumed to be 0.05 (a relatively high 

figure is chosen to improve the visual impact of the root locus plot) . 

From Eq. (298) of Reference 1, the poles and zeros of the open-loop 
transfer function G^Cs) may be obtained. With normalized to unity, 
these roots become 

?! 2 = 0,0 ^ 12 ^ ^ ^ 

P3 = 0.111 ± i 1.49 

For Figure 3 the controller (feedback) transfer function 113(5) has a 
pair of complex poles, a real zero and a real pole of larger magnitude 
(so it might be called a lead-lag system "with complex poles). The influence 
of flexibility on the locus is not surprising, but the "crossover gain" 
which marks the transition to instability is drastically reduced by the 
flexible appendage. Whereas a rigid vehicle would become unstable when 
the gain reaches 0.13, Figure 3 shows instability to result when the gain 
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Figure 3. Test vehicle root locus for lead-lag control with remote complex poles 
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Figure 4, Test vehicle root locus for lead-lag control with near complex poles 
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Figure 5. Test vehicle root locus for lead-lag- lag control with complex poles 
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Figure 6. Two-mode flesdble test vehicle root locus for lead-l'ag-lag control 



becomes 0 . 05 . This is a reduction for the flexible vehicle to 38 % of 
the cross-over gain for the rigid vehicle. 

parenthetically, it may be noted that the gain K employed in these 
root locus plots is related to the gain which appears as a factor in 
the feedback transfer function by a factor which is uninfluenced by the 
flexibility of a single mode. Thus changes in root locus plot gain K 
due to flexibility can be interpreted directly as changes in the actual 
feedback controller gain . 

Figure 4 is, like Figure 3, based on a lead-lag control system with 
a pair of complex poles. The poles of Figure 4 are somewhat 

closer to the original than in the previous example. The cross-over gain 
for the flexible vehicle appears to be 20^ higher than for the rigid 
vehicle, whereas this gain is 62 ^ lower when the complex pole of H^(s) 
matches the structural frequency, as in Figure 3. Thus generalizations 
regarding the influence of appendage flexibility on control system 
stability would appear to be dangerous, and superficial physical interpreta- 
tions unwise. 

Figure 5 portrays the rigid vehicle response to a "le ad-lag-lag " con- 
troller, with a negative real zero and tv;c more remote negative poles (one 
of which is off the plot at -2.0). The cross-over gain is K = 0.31. 

Figure 5 illustrates the influence of the first mode flexibility of the 
appendage, which reduces the cross-over gain to 0.13, 42^ of its original 
value. Figure 6 indicates the influence of the first and second modes of 
appendage vibration in combination. 

For a "pure gain" direct feedback control, root-locus plots (not shown) 
can be deduced by inspection. It can be established by use of the Routhian 
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array that such control systems cannot be destabilized by appendage 
flexibility. (See Ref. 1, pp. 72-74 for details.) 

For this idealized test vehicle there were generated 35 root locus 
plots, of which those preceding are typical. With a standardized program 
(written by E. H. Kopf and R. Mankovitz), each plot required only 10 to 
4o sec of IBM 7094 tine, depending on the number of plots requested in 
a given run. For the single-axis linear response of a vehicle as simple 
as this test vehicle, the conclusions are quantitatively precise. For a 
vehicle of realistic complexity (as treated in subsequent appendices), root 
locus plots nay lose quantitative validity, while preserving the qualitative 
significance normally sought in preliminary design studies. 

APPENDIX B — THE TOPS AUTOPILOT ; 

The Thermoelectric Outer Planet Spacecraft (TOPS) configuration (Figure 
7 ) is dominated by a(4.27ra) dian parabolic communications antenna and a 
bank of radioisotope thermoelectric generators (RTGs). In addition, a 
pair of telescoping booms are required to provide separation for sensitive 
instruments (magnetometers, plasma wave detector, etc.). Much of the 
remaining scientific instrumentation has been provided viewing area (around 
the antenna) and separation from RTG radiation by mounting it on a large 
foldout structure opposite the RTG foldout boom. Finally, central to the 
craft is an electronic equipment compartment carrying the autopilot's 
attitude sensors (high-gain gyros) and actuator (gimbaled engine). It is 
this portion of the vehicle that will be considered the rigid body to which 
the flexible appendages (RTG, science, and magnetometer boom.s; antenna) 
are attached. 
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Figure 7. Thermoelectric outer planet spacecraft 


While the hybrid coordinate formtiLation is also the ideal approach 
for the TOPS cruise attitude control system because of the use of momenturti 
■wheels, the discussion here is limited to initial investigations of the 
autopilot 'Which maintains vehicle attitude during the trajectory correction 
motor’s thrusting periods. The gimbaled engine provides attitude control 
about the pitch and yaw axes, while roll con-trol is maintained by the 
roll gas jets, which, in the cruise control mode, are normally used to 
desaturate the roll wheels, pitch, yaw, and roll momentum wheels are 
switched out during the trajectory correction phase along with their respective 
cruise optical sensors. Rate and position sensing is provided in all 
three axes by high-gain, rate-integrating gyros. 

Preliminary efforts have attempted to determine the feasibility of 
using an autopilot control loop originally designed for a Mars orbital 
vehicle on TOIS as well. A linearized and simplified version of the pro- 
posed autopilot yaw axis (or pitch) loop is sho'wn in Figure 8. Among the 
assumptions made to develop the simple, single-axis system model was that 
gimbaled engine reaction torques on the craft would not be considered, i.e., 
the engine is assumed massless. Also, as discussed previously, the influence 
on 9 (pitch or yaw) of constant thrust along the vehicle roll axis is ignored; 
only the effects of the torque applied by the gimbaled motor and the resulting 
rigid-flexible body interactions are included. 

Of course, the three-stage design process to be executed for the TOPS 
autopilot system must be based on the existence of a detailed hybrid coordi- ' 
nate structural model of the craft pictured in Figure 7, Such a model was 

developed relatively quickly (2 weeks) using 20 discrete sub-bodies to 
* 

Thanks must go to Messrs. Robert Bamford and Craig Heiberg of Division 35 
who modeled the structure and performed the necessary computer calculations. 
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Figure 8. Linearized and simplified yaw autopilot loop 





approximate the flexible members. The result of computer calculations to 

transform the discretized structural model into hybrid coordinate modal 

data is summarized, in part, in Table 2. Shown are the first ten modal 

frequencies obtained and corresponding elements of the matrices 6 and 
T 

0 Approximately 30 modes were available in all from the computer 

generated data. 

It is at this stage that one must begin to weigh the importance of 

the modes and, tentatively at least, truncate the modal data. The first 

five modes, grouped near 1 Hz may be traced to the antenna and magnetometer 

booms. One can see from the size of coefficients in 6 that pitch and yaw 

axes are predominantly affected by modes' 1 and 4 while roll largely sees 

modes 2, 3, and 5. Based on very simple modeling of the craft as a rigid 

body hinged to two rigid booms, natural frequencies of vibration above 

4 Hz were shown to have little effect on autopilot stability. Therefore 

modes beyond 5 'were removed from the dynamic model for all phases of the 

analysis including the detailed computer simulation. An additional excuse 

for deleting the significantly higher frequency modes is the dramatic 

improvement in numerical integration speed if digital simulation is employed. 

Equation (298) of Reference 1 can now be of use in examining the 

linearized, single-axis control loop-. For both pitch and yaw, modes 1 and 

4 appear to dominate in the 5 matrix, so that 9 and 0 could be approximated 

y X 

by equations of the form 


8y( = ) 


[I s" 

yy 


— y y.. 


(s^ + + cr^) 


- s 


6 6 - 

y y 




■] 
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Table 2. TOPS hybrid-coordinate structural model data 


Mode 

i 

^i« 

Hz 

T 

5 MZgQ (kg-m) 

6(kg-m^) 

X 

y 

z 

X 

y 

z 

1 

0.74 

0.2839 E-3 

0.3557 E-3 

0.1820 EOl 

-0.9479 B 1 

-0.1559 E02 

0.5230 E-2 

2 

0.75 

0.1392 EOl 

0.1565 EOl 

-0.6839 E-3 

-0.5334 EOO 

0.4858 EOO 

0.1802 E02 

3 

0.76 

0.1319 EOl 

-0.1610 EOl 

-0.4712 E-4 

0.5523 EOO 

0.4553 EQO 

0.1700 E02 

4 

0.76 

-0.1766 E-3 

-0.2151 E-3 

0.2348 EOl 

-0.1216 E02 

0.1172 E02 

-0.1953 E-3 

5 

1.16 

0.1235 EOO 

0.5553 E-2 

0.6532 E-2 

-0.2891 E-1 

0.1992 E-1 

0.6242 EOl 

6 

3.85 

-0.1216 E02 

-0.13365 EOO 

-0.6915 E-2 

0.2269 EOO 

0.8294 EOO 

-0.3575 E02 

7 

5.02 

0.1782 EOl 

0.1965 EOO 

0.1079 EOO 

-0,8941 EOO 

0.5455 EOl 

0.1501 EOl 

8 

5.66 

-0.2009 E-1 

-0.1075 E-1 

-0.8343 EOO 

0.1163 EOl 

0.2636 EOl 

-0..9897 E-1 

9 

5.66 

0.9229 EOO 

0.7903 EOO 

-0.3320 E-1 

-0.1688 EOO 

0.3133 EOO 

0.3625 EOl 

10 



5.69 

-0.7921 E-1 

0.4353 E-1 

0.1075 EOl 

-0.1492 EOl 

0.2003 EOl 

-0.2894 EOO 


I = 


1544.9 

- 2.29 

- 2.76 


- 2.29 - 2.76 

471.9 - 34.9 

- 34.9 1714.4 


kg - 562.6 kg 



(The equation for 0 ^(s)/t^(s) is identical in form). 

With the substitutions I , I , 8^, 8 8^ =(1544. 9» 471.9? -15.6, 

XX yy y -X. y 

-9.479,11.722, and - 12.16 Kg-M % respectively, and = ^.65 and 

4-. 76 rad/sec 5 respectively, and = 0.005, the resulting transfer 

functions for yaw and pitch are 

(s^ + 0.0465s + 21.62) (s^ + 0 .o476s + 22.8) 
1139 s\o. 405 s + o.o66s^ + 31. Is^ + 2.09s + 493] 

2 2 
(s 4 - o.o465s + 21.62) (s + 0.0476s + 22.8) 

1139 s^[ 0.887 s^ + 0.089 s^ + 4l.9s^ + 2.09s + 493] 

The open-loop poles and zeros of the linearized, single-axis autopilot 

system of Figure 7 may now be plotted as shown in Figures 9 and 10. The 

dashed lines on Figures 9 and 10 indicate the closed-loop root locations 

under the ideal, totally rigid spacecraft condition where the control 

parameters are nominally specified as; Xn ~ 1-77 sec, = 0,111 sec, 

Cq - 0*35, Kp = 2.2, oj^ = 88.0 rad/sec, = 3.47, Tp = 20.0 sec, and 

(u = 138.2 rad/sec3 gimbal servo poles; -938.66, -20.35; gyro poles; 
s 

-30.8 ±82,4 i; autopilot poles: -9.12, +0.059; gyro zero: -O.565; and 

path guidance zero: -O.05. 

Loop gain (DC or Bode gain) at the point of marginal stability is 
10.0, On the other hand, in Figure 9 the yaw axis loop with approximated 
flexible spacecraft dynamics shows a drastically changed root locus with 
a critical gain reduced to 2.9, While the pitch-axis root locus for the 


Ty(3) 


0 (s) 
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flexible case also departs significantly from the rigid case, critical 
gain reduction is only to 8.5. The fact that the total pitch axis moment 
of inertia is about 3.3 times as great as yaw's obviously was responsible 
for the difference in location of the flexible appendage pole -zero contri- 
butions to each axis and the effect on critical loop gain values. 

An examination of the rigid-body autopilot system reveals closed-loop 
bandwidths of about 0.3 and 0.8 Hz corresponding to DC loop gains of 
1.0 (0 dB) and 3.16 (10 dB), respectively. It is clear that the first 
four appendage modes have already entered the control bandwidth at the 
latter gain level and, in the case of yaw, caused instability. 

One might expect that, if additional modes are used in the initial 
single-axis root locus analyses, critical gain values would decrease even 
more. Further, the fact that interaxis coupling might require even more 
substantial loop gain reductions to maintain stability suggests an 
eigenvalue analysis of the complete system. The TOPS autopilot may be 
arranged into state variable form to give a 30 x 30 system state matrix 
including pitch, yaw, and roll angles, five flexible appendage modes, and 
the various control dynamics. Figure 11 shows that when the detailed, 
coupled case is considered, yaw loop stability deteriorates further and 
essentially controls the stability of the entire system. 

A digital simulation was programmed and included nonlinearities due 
to (1) gimbal servo drive amplifier saturation, (2) the saturation 
characteristic built into the compensation block to prevent reaching 
mechanical gimbal stops, and (3) the roll axis bang-bang control loop. 

(it did not include gimbal actuator stiction.) The level of detail was 
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Figure 9. TOPS autopilot yaw axis root locus 
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Figure 10. TOPS autopilot pitch axis root locus 
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such as to include effects of gimbal mounting errors and center- of -mass 

displacements due to appendage deformations, and the hybrid coordinate 

formulation represented by Eqs, (224) and (2l6) of Reference 1 was 

implemented in detail, yaw-axis response is shown for three cases in 

Figures 12 and 13 to illustrate the computer simulation results. In each 

case, the autopilot controller and spacecraft are responding to the engine 

turn-on transient ( 445 H ) under the conditions of an initial pitch 

axis center-of-mass offset of (.05 cn) Five flexible modes are used in 

each case in Figure 13; for gains of 2.2 the flexible vehicle proves 

to be unstable (Figure 13a) as the eigenvalue analysis predicted. Reduction 

of K to 1.0 results in a stable response (Figure 13b), which also 

x,y 

agrees with eigenvalue predictions, 

APREKDIX C — THE ATTITUDE CONTROL OF A FLEXIBLE SOLAR ELECTRIC SPACECRAFT : 
Introduction ; 

During the cruise-low thrust phase of the mission of the solar electric 
spacecraft illustrated in Figure l4, one means of providing three axis at- 
titude control is the combination of translation and gimballing of the 
electric propulsion engines. This concept will be analyzed in the sequel. 
Inclusion of the elastic effects of the flexible solar arrays was ac- 
complished by utilization of the hybrid coordinate procedures described 
in Reference 1. 

Figure l4 shows the spacecraft to be analyzed. It is composed of a 
rigid central body to which are attached two flexible roll out solar 
arrays A, and an engine cluster E. Unit vectors bj^, ^ and b^ are fixed 
in Figure 15 shows E in more detail. The engine cluster consists 
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of a large plate to which are attached three electric propulsion engines 
and a unit containing translating mechanisms. Two of the engines, and 
can rotate about the diagonal passing through the ~ ^ quadrant. 

The plate and the remaining engine are denoted B 2 . E can translate 
parallel to the direction and B 2 can translate parallel to the b^ 
direction. Figure 15h shows the cluster in an arbitrary configuration, 
point 0 is the mass center of the entire spacecraft for the case of no 
engine translation or gimballing. 

During the cruise low- thrust phase, when E translates parallel to 
a control moment parallel to bg is exerted on the spacecraft. Translation 
of Bg parallel to b^, will produce a moment parallel to Third axis 

control is provided by the gimballing of engines and Bj^.. This generates 
a moment parallel to 

Figure I 6 shows E in even greater detail. In addition to Bp, B^j and 
Bjj^, bodies B^, B^ and B^ are shown. These are associated with the mechanisms 
which translate E. and Bg are drums with axes fixed in the translating 
mechianism unit. Two tapes are wrapped about Each has one end attached 

to B^. The remaining ends are fixed in Bj^. Hence, when rotates, one 
tape winds up on the drum while the other unwinds. Consequently, E 
translates along the shafts parallel to b^, A second assembly consisting 
of Bg and two tapes accounts for translation parallel to except that 
the tapes are attached to Bg, instead of B^. is the unit containing 
Bj-, and the motors which drive them. 

5 6 

Differential equations of motion for the entire vehicle, the appendages 
A, and the rigid bodies Bg, will be stated. 
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Vehicle Equations; 


For the assembly of bodies shown in Figures l4 and l6 the vehicle 
equation is, from Reference 1, 

a ot> 0, 

T = + [c+2ujxc 




(OX (wxc) + tOX£] XC + 



I £ 


o 

E 


dm 


( 1 ) 


where 

T - torque applied to vehicle. 

I - vehicle inertia dyadic. 

m - angular velocity of in inertial space. 

- total vehicle mass. 

£ - mass center shift vector. 

J - an inertial reference frame. 

g - position vector of a generic point of the vehicle with respect 
to 0. 

dm - differential element of mass. 

(*) - derivative of a vector or dyadic in b, or the derivative of a 
scalar. 

(°) - derivative of a vector or dyadic in b, a reference frame in which 
the unit vectors b^, b^, b^ are fixed. 

In more detail, the matrix form of the vehicle equation is 


T = [I 


* 


i=2 


+ T 


IT. 


lV)]cB 
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+ U)[I 


+ T. (-M.p p + T I T )]w 
i=2 ^ 


+ i;2M(EgQR + p)VRq^M SgQ- - p'^ (Mqf ^ 



t 

P 


rgo(qM - Mq)S^^}u> 


OE" 


c- 


1=^2 




itvf'] 


i=2 


oj + TTL(c +2 5^6 + 03 + (ac) c 


+ 


7 


s Cm (?? 

i=2 


i*..i* 

r 


+ ar p 


i*.i* 

p 


) 



X 






T I T (u } 


(cKR) 




+ R t^^Mq + R ® 7^QMq 




+ T' 
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- ’’’ ® ^0E'“ 

- (2 T:^g>w+ 'iT{2[M(-gQR + p)]^q 

- Rq'^M 

+ {2[m(Sj5qR + p)]^q - Rq'^>i r^o - 

- p‘^(Mqf^ - CMq)"^ P^'^h (2) 

where 

T - applied torque matrix, 

'J'T 

I - inertia matrix for and A in the undefomied state. 

M. - mass of B. . 

1 1 

-x- 

p - position vector of the mass center of with respect xo 0, 

- transformation matrix from principal axes of to principal 

axes of 

- inertia matrix for 

u) - angular velocity of B-j^ in . 

M - the nxn mass matrix for A where n is the number of sub-bodies 
comprising the model which approximates A. 

R - position vector of a point Q on the appendage - interface 
with respect to 0. 
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0 - position vector of the mass center of a sub-body of A with 


respect to Q. 

q - deformation coordinate vector for A. 

SgQ - the matrix [EOEO. . ,E0] of size 6nx3 where E is the 3x3 identity 
matrix. 

T 

j; _ - the matrix [0E...OE] . 

'OE 

c - mass center shift vector. 

T 

~ - when applied to a 3 x 1 vector, v = C‘V 2 _'V 2 V 2 ] , v means 

To -V3 vj 


T 3 0 -v^ 


^ 2 n i 

when applied to the 3 ^ 6n matrix v = [v Ov 0. . .v 0] , where v^ is a 3 x 1 
vector, V means 
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T 

( ) - the transpose of a matrix. 

( - if V = [v^Ov^O. , .v^O]^, then 

= [v^Ov^O. . ,v”o] . 

An expression for c is 


c 



( 3 ) 


where 6^ is the displacement vector for 

If engines and are thrusting at the constant magnitude F, 

and there are no disturbance torques exerted on the vehicle other than from 
the thrust, then 


T - m)r^ + ('m - m)r^ 


rn ^ 




+ (7H - m)r^ - E^QMq]} ^ 


FS3 

-FS3 

/2FCj^ 


+ [p^*^ + ~ [( Tn *- m)r^ + ( TlX - 


- 4Mq]r^ 


Fs, 


-Fs, 


|/2FCi 
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where 


i^O ^ 

p - position vector of P^, the mass center of with respect 

to 0 for the nominal state. 

7 

m - y. M. . 

i-2 ^ 


m - M2 + 

- [-^or^ 0 0]*^. 

r^ - [0 -rtig 0]^. 

Ck'j^ - the angle through which turns relative to B^. 

s ^ — s In > i“2 ^ ^ 7 • 

c - cos O’ , i=2,...,7. 

1 3 j 1 

r - the rccJTUs of ci Bg. 

Appendage Equations ; 

The appendage equation is taken directly from Eqs. (9^) and (95) of 
Reference 1. 

* t2MC(Eg^r - 1 


+ MteQgtu) + (Sq^O)) M - + D'}q 


+ {M(rOE<i) - (%eW) - (SQjfOrCMSQgw) 
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* ^ * n']q 

= - MEol” - (Soi® 

- (R + ‘3?)uj + e - ^ w + aJ ©(R + e)] 

- p- Ego* - (r^o”) 

Where the nei; symbols mean 

- damping matrix. 

- stiffness matrix. 

qi 

ex - acceleration of 0^ in J expressed in terms of {b] = 2— 3^ 

X - forces and torques, other than from, structural interactions v?hich 
are applied to the sub-bodies. 

e - the mass center shift vector associated with i=2j...,7. This 
is the second terra in the expression in Eq. (3) for c. 

Rigid Body (i=2, ...7) Equations; 

The units . .B,^ constitute the translator. Equations for the co- 
ordinates associated with its translation and the giniballing of B^ and 
are required. It appears as if use of Lagrange's equations provides the 
most convenient means of determining these equations. The standard fom 
for these equations is 
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( 6 ) 


(it ^9^-' ~ ay^ 


% 


where L is the Lagrangian for the spacecraft, y^, is a generalized coordinate, 
and is the generalized force corresponding to y^^. L is 

L = K - P (7) 


where K and P are the kinetic and potential energies, respectively, of the 
spacecraft. Equations (2) and (3) govern the behavior of u) and q. In 
using Eq. (6), then, ve \?ill restrict attention to coordinates other than 
oj and q. 

The following assumptions are made for xhe translator assembly (refer 
to Figures I5 and 16) : 

1. For the translator assemblies, shafts connecting either or 
Bg to the motors driving them have non-infinite stiffness. 

2. Won- infinite stiffness is associated only with the drive shafts 
which cause and to girabal. 

3. All drive shaft bearings are not smooth. 

4. The mass centers of and lie on their gimbal axes. 

5. The translation shaft bearings are not smooth. 

The nonlinear equations of motion for and are 

4^3 ■^75“ (“1 * “2^^ ■ ^^3 ■ 4)^75“ ■ “2^ 


- sgo, 




23(0,5 


“2) =3“3^- 


- Qf^) = 
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- %K - ^2^ + cj^to33 *- y.^{a^2 ~ "4) -V4 (9) 

(mr^ + - mrE^^GX - fiiu)^ + I^di^ 

- mr + 103^^) - (tu^ + (x 4 )mr^a^ 


- 2 - 2 . 

+ mr “^3*^6 


3T 2t 
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/A A \ / A A \ 

+ ( tu ^ ^ ^ 2 ^X 


^"“ 5"3 ■" t^l "3 ' "irT"' °'^l "2 * bT ’'^“' 5'"2 


2m 2 2 


— 2 2 2 
mm2 m .A . m .A m r, ^A J^^, 

-r OoOi-, - tt— rc^oj^ + - -:c;^Coj 2 C;j^ - rUj_C 2 Ju )3 


— r 0 ) 20 ^ 3 ^ - - W 
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2 . 


( 10 ) 


/- 2 ^ 6 x.. “ „ 2 T^. - . , ^ 6 . 

(mr + ^3 ^^6 ” ^ ^ ^ 3'^3 


+ mr ^{ tb ^ a ^ + (u^a^) - (cu^ + + mr'^tio^ujgGr^ 


2 . 2,-2 - 2 


— It 3t • 2 • 

- hhu 2 (U 3 - mr(E tu^ - E ''" *^3*^5 
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— r (c .^3 f c.^3) + :j5i;-r(«)3C3 + .-.303^ - 0,3=3 


48 


JPL, Technical Memorandum 33-493 



where 


.A , m mr/ \ 


- mm3(«)2‘=3 - ^3^^ ■" 


2m^ 2 2 . 2 2 


mm2 , nur , — _ _ ... ^ _ 

— r + ^2.1.3 + :j^r 0.^3 + 


iri El , A m m .A m mr, A Av 

?a ~ m m ' 3 2^ 


— 2 ~2 

mmr/A Ax 2mnr 2 2m.m 2 

m ^'■“l^2 “ ‘^2^1 ^*"1 ” “ 'w~^^l 


I'F/ N . 2. 

" " 7T^^3 " V ""6 


( 11 ) 


It - moment of inertia of B, for the axis parallel to b 
J 1 ^ -j 

o. - b. • (Uj i = 1,2,3. 

1 -— 1 — 


- spring constant associated with the translating or gimbaling 
mechanisms. 

- damping constant associated with the translating or gimbaling 
me chanisms. 

eIT^ - [1 0 0], [01 O], [001]. 

c^ - measure number of mass center shift due to flexibility. 

i* i*, , i*0 ^ i* 

li - where is the component of p , the value of p 

for the nominal state . 

- torque applied to B^. 
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a . - gimbal drive shaft input angle, 

Gi 

~ translating drive shaft input angle. 

Linearized Equations of Motion ; 

It can be shown by substitution that 

0, = o'i = X = cv^l = = ‘^Tl = ®T2 = ° 

X = 0 = Oq (12) 

q = qo = - 

is a solution to Eqs. (2), (5), (8}-(U). Determination of the linearized 
variational equations associated with this solution and introduction of the 
active control equations establishes the base for a stability analysis by 
means of root locus procedures or eigenvalue analyses. 

The variational equations are obtained by substitution of 
CO = to' 

i i/ 
a = a 

*^Gi^ ^Gi’ “Ti ^ °^Ti 

X = X + x' 

o 

0 = 0 + 0 ' 

o 
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into Eqs. (2), (5), (8)- (11) and linearizing in the primed variables. 
\ is set equal to zero. 

The appendage equation takes the form 

U'-i' * B'4' + K'q' = t 


7 


•where 






In all of the following discussion, the primes -vyill be omitted from the 
linearized system variables . If the modal coordinates Ti are defire d by 
q = j^T], -t/here f is the matrix of eigenvectors for 


M'q + K'q = 0 

and if ^ and T| denote N x JI truncated versions of ^ and T| "where K < n, and 
if 6 and A (hereafter referred to as the rigid-elastic coupling terras) are 
defined by 

6 = - (Sqj. - 


= ^mSbo 


then the appendage equation becomes 


• • • _o 

E71 + CT Tl + a f 


6e + ^ (oj ^ 
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■where E is the N x N unit matrix, a is the N x N diagonal matrix of eigen- 

•values or mode frequencies, Q is the N x H diagonal matrix of damping ratios, 

2 12 
F is the thrust vector for the center engine, and r and r are vectors 

1 2 

describing the b and b migrations of the translator from its nominal 
position 


r^ = C-rcv^, 0, O]"^, r^ = [0, -ra^, O]' 


It "will be assumed that = ”^3’ appendage equation becomes 

(11;) 


ET| + 2(; CT 71 + CT^ = S*9 + + mr^) 


The torque T exerted on the vehicle may be expressed as 

^ ~i i 
■p = ^ s F 

1=2 

■where F^ is the thrust force supplied by i=2, 3, 4, and is the position 
vector of the point of application of F^ with respect to the vehicle mass 
center. It can be shown that a more explicit non-linear expression for 
T is 


^ r 2*0 m- m 1 - m - m 2 1 T 

T = IP + r + ^r v Mql 


Trt 


in EO 


0 

0 

F 




7U 


m. 


1 j 1 


Fs 


3- 


-F£ 


/2Fc, 


5Z 
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IflL- 


m 


m 1 
— r 




m- m 2 

m 


TTl 


1 

7T 


-Fs^ 

/2FCi^ 


Linearization of this expression and of the right hand side of Eq. (2) 
about the solution of Sqs, (12) produces the following linearized vehicle 
equation 


10 + Kr + icr + 


, -, 6 .. 6 

+ I /y 




4*0 ~ 

P ) 



- - m)r^ + (TYl - m)r^ - A^] 


(15) 


New symbols are defined as follovfs 

I - vehicle inertia dyadic for the nominal state 

i*0 i*0 

P - E M-P 

i=2 


v>l p' 


0*0 


0 = 2 , 3,4 


Truncation of all but the first six modes ■57ill be adhered to in the 
sequel. Thus 5 and A are 6x3 matrices. A structural analysis of the 
entire vehicle reveals that only three of the elements of 6 and two of the 
elements of A are non-zero. In particular, 6 and A have the forms 


~0 

6 12 

0 

5 A = 

0 

0 

0 

0 

0 

^23 


0 

0 

0 

0 

0 

0 


0 

0 

^33 

0 

0 

0 


0 

^^42 

0 

'51 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 


JPL. Technical Memorandum 33-493 


53 



1 2 

Substitution of these expressions and of those for r and r into Eq. (l4) 

/ 

produces the following six scalar equations of motion for the modal 
coordinates 

\ ^^1 ^ 1 % "" ^ 12®2 

^2 ^^^ 2^2 '^ 2*^2 ^ ^ 23®3 

(16) 

** *2 

’I ■" “ ° 

where < 7 ^ and are the frequency and damping ratio respectively for the 
ith mode. 

Figure 17 shows the mode shapes and lists the frequencies corresponding 
* 

to the first six modes. 

The vehicle matrix equation (I5) is then equivalent to the three scalar 
equations 

Provided by Dr. E. Weiner of Division 35 > "who modeled the structure and 
made the necessary computer calculations. 


'Z 
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^ 0-1 = 0.02129 hertz 



cr^ = 0.02376 hertz 


^ MODE 2 
7 0*2 = 0.02141 hertz 



MODE 5 

o'5 = 0. 02838 hertz 


MODE 3 

0-3 = 0.023589 hertz 



MODE 6 

cr^ “ 0.028315 hertz 


Figure 17. Solar Array Mode Shapes 





- 'l2% = 3^^rF«5 


^3 ^“5 ~ " 4 A 


- 623T12 = ^2p^%a^ (I9) 

where is the moment of inertia of the vehicle for b^ and for the nominal 
state, and are the measure numbers of K and K, and F is the 
magnitude of the thrust force exerted by. any of the three electric propulsion 
engines . 

The assumption that = -a^ implies that the equation of motion for 
may be ignored. The linearized forms of the remainder of Eqs. (8) -(11) 


^^“3 ■" * S“3‘ '"3^1' "3^ “ ° 


1 2 2 

/-t 5 ^ 2 , m^ 2 ^ r s.. 2^ . 


- - 45, - 0 


2 ^3 3 


6 ~ 2 mm^ 2 2v.. 2, . ^ , 
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2 — 

m m^2.. 

m2 "5 


(f^r - 


m 


- - iTr)'e3_ + - 0 (22) 

Control : 

The linear approximation to the control system block diagram for each 
axis is shovn in Figure I8. A celestial sensor detects the value of 0^, 

The transfer function for the sensor is Kgj^/(Tgj^s+l) , A constant 
is the transfer function for a voltage controlled oscillator while the linear 
approximation to a stepper motor has the transfer function • The 

remaining part of the control system is the compensation network. The 
output from the stepper motor is the desired translation ^^2^ 

or the desired gimbal angle of , The sign associated with V^. is determined 
by the sign associated ;;ith the translating or gimballing controlling term. 

A block diagram of the dynamics and control of the entire vehicle 
is shown in Figure 19. Coupling between the three axes is indicated. 

Differential equations associated with the control system are obtainable 
from Figure I9. They are 


1 = w 2. 3 


(23) 




^T2 ^42*S/2^ ^S2"\2^ 




(2lt) 
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U1 

00 



TRANSLATION ra^, 

OR GIMBAL ANGLE a g 


Kli Km, 
Tli S + 1 



l_ 1 l_ 




li>p] 








































V, 




FI 




{ Kqi®i + (l+K, 

SI 


'21 


)^c 


'^21 Tgi' SI 


Tii+ T2i 

^^21 ^ IViK^i^V FI 


- (1 + K„. + 


)v^i3 


21 FI 


'F2 




'^12'^ 22 


'22 

_ '*''<52 

^S2 


K 


+ (Too + 


'’■ 12 '^ '*’22 


)v. 


V ( 25 ) 


22 ^2^12^12^^2^ 1’2 


+ (1 + K„ + 


)V^o3 


V. 


22 IV2^12^I2^2 ^2 

-^3^13^^3r'^13 


- 9 + (1 + K,,- ;;S)v 

F3 ^^13^23 ^S3 ^ ^S3 


** (Toq '*' 


"’ 13 " -)v. 


23 Wi3"m3^3 '’3 


- (1 + K_ + 




23 IStjKls'Sb'Sq « 


stability Analysis ; 

If the coupling and translator — gimbal dynamics are ignored in Figure I 9 , 
then the block diagrams shown in Figures 20-22 result. The drawings at 
the bottom of each figure indicate the mode or modes which are excited by 
rotation of the spacecraft about the axis to which the block diagram cor- 
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responds. As indicated by Eqs. (l6), the third and sixth modes of 
structural vibration decouple from the remaining equations and consequently 
have no effect on attitude motion, at least in the linear sense. Uncoupling 
of the equations transforms the attitude control problem into one that is 
amenable to root locus analysis. 

The open loop transfer functions for the three axes are 

' ■- 


+ N^s + W^)(s + 1 /t 2 i)(s + " 


2 34 3 2 

6^l)TgiS (s + D^s + DgS + Dj^s + Dq)(s + 


VTsi)t=^ ^ ^ 


+ + 


^ii^aNrAi 


'^21'^u ^2i’’‘n 


(1 + K23_)3 


(26) 


^2^2 




^^2”'^ 12 ^'^32 


ic^m 


2 s CT^)(s + 1/Tp,)(s + 1 /t,^)/s3(s2 + 


<22''' f '12' 


21 


I CT^ 

; — ^ Vt32){s^ + 

±2 - I2 - O12 
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+ _1_ + ^ 


T 22 TJ 2 '^ 22^12 


(27) 


' ' (^ - «23>^S3 


lA,3)(s.lAy)/s3(s=.!^s.!fV) 

^3~®23 ^3'^ 23 


(s + 1 /Tg 3 )[s + ^^^ 3 ^ 13^3 


. ^ H- SsSe^a . ic^)) 

^23^13 ^23^13 


(28) 


where is the for\?ard loop transfer function and is the feedback loop 

transfer function. Expressions for K-, i=0,...,5 and for D , j=0,...,3 

1 J 

are 


" 2(^2^0'2j^ + 


h 


2 2 


* ‘<hCfh°5 


+ 3('rtX-m)F 
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N- 


= + Qfl^) + 

3 


6FmA^2C^5 


H. 


' ”l’’5 * 


iT3m 


— 2 2 
3FmA^2'=^^ 


6(1U - m)FcTjO 




"o - 


3(m - WFg^cr^ 
^3^ 


(29) 


and 
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D-. = 


Vc2(cu<>4 




V »51 


D, 


+ <.5 + ‘*C4C^4C5 


V =51 


.2 2 

Sl“■l(^ 




”1 “ :;^ — 12- <^if“5 * 4) 


^1 " ^51 


®0 = 


T 2 2 

I^c•, a^ 

145 

V=51 


(30) 


The values selected for the parameters in these expressions as well 
as for the parameters in Ejs. (26)- ( 28 ) are 


^li 

= (328 volt seconds/meter) 

i 

= 1, 

2 , 3 

^2i 

= 4.5 

i 

= 1, 

2 , 3 

^i 

= 80 

i 

= 1, 

2 , 3 

% 

= K^2 = ((^.7 X 10 ^meter/pulse 




"Sb 

= 10 ^ radian/pulse 

i 

= 1, 

2, 3 


= 4 o pulses/volt second 

i 

- 1, 

2 , 3 

^23 

= il 28.2 Kg ^ 2 ) 




^12 

= (128.3 Kg m^) 
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^51 = (9.609 Kg m^) 

^42 = "(13.69 Kg m) 

"i = (.006 Kg'^ m) i = 1, 2, 4, 5 

= (.02479 X 2n Kg^ m/second) 

^2 = (.02493 X 2tt Kg^ m/second) 

°4 = (.02767 X 2tt Kg^ m/second) 

_ >- 
5 = (.03305 X 2"t Kg^ m/second) 

T 3 

li = 10 seconds 

T 2 

2l = 5 X 10 seconds 

'Si = 1 second 

F = (.09 Newton) 

= (442.7 Kg m^) 

- (12777.0 Kg m^) 

= (12408.0 Kg m^) 

= (600.0 Kg) 

m = (54.41 Kg) 

III .= (45.35 Kg) 

= (-41.47 Kg tn) 

= (-34.55 Kg m) 

The sensor gains K . have not been specified and will serve as variable 

o X 

parameters ypon which the attitude stability depends. 

The root locus for the b^ axis is shown in Figure 23. Open loop 
poles and zeros corresponding to the b^ axis case are listed. The poles and 
zeros having non-zero imaginary parts are associated with the fourth and fifth 
flexibility modes. A real zero at 0.084324 is due to the so-called "tail- 
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- 0.2 

■ LATJV 

0.0 0.2 

0.4 

0.6 0.8 


REAL 


Figure 23, bj Axis Root Locus, Fourth and Fifth Mode Flexibility 


vO 



TB-gs-clDC effect" or the term *<^2^6 inertia torque 

exerted on the vehicle when the engine cluster unit translates in the 
direction. Such translation activates fourth mode oscillations. Figure 
2k is an enlarged version of Figure 23 and shows detail near the origin. 

For very small gain, i.e., very small "the locus emanating from the 

multiple pole at zero initially passes into the right half-plane. Thus, 
instability is implied. As is increased, the locus crosses back into 
the left half -plane and finally returns to the right half-plane for suf- 
ficiently large Kg^. That part of the locus vjhich corresponds to flexibility 
traces a path from the left half-plane into the right and then back to the 
left half -plane. The range of for which all parts of the locus lie 

in the left half-plane is 0.50 < < 23.0. For the case when flexibility 

is ignored and the vehicle is assumed to be rigid, the stability range for 
Kg^ is 0.52 < Kg^ < 23.52. 

Figure 25 shovjs the root locus for the axis parallel to bg. As before, 
the poles and zeros .not on the real axes are flexibility poles and zeros . 
Also, for very small gain, the locus traced from the multiple pole at zero 
passes into the right half -plane. As Kq 2 increased, the locus moves 
into the left hadf -plane, then to the right and finally to the left of the 
imaginary axis. A part of the locus associated with the flexibility pole 
begins in the left half-plane and soon moves to the right half -plane. 
Stability for the axis is guaranteed for I5 < < 201. For comparison. 

Figure 26 shows a root locus for the b^ axis when the vehicle is assumed 
to be rigid. The character of the locus in Figure 26 differs considerably 
from that in Figure 27'.' It was determined for the rigid model that the 
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Figure 2h, b, Axis Root Locus 
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Figure 25 . Axis Root Locus, First Mode Flexibility 
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Figure 26. Axis Root Locus, Rigid Vehicle 
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Figure 27. ^3 Axis Root Locus, Second. Mode Flexibility 




range of vtiich vill ensure stability is l 6 < 276. Thus, the 

flexible nature of the solar array is detrimental to the attitude stability 
of the bg axis. 

The final root locus corresponds to the axis parallel to b^. There 
is no "tail-wags-dog" term for this axis, and this is reflected in the 
absence of a zero on the positive real axis. As in the case of the bg 
axis, there is only one flexibility pole and one flexibility zero. Values 
of Kg 2 "Which stable attitude motions are guaranteed are 11 < < 221 

for both the flexible spacecraft model and for the rigid model. 

Time histories of the attitude angles 0^, 0g, and 0^ in Figures 

28-33. The CSSL simulation language, developed through the support of the 
Guidance and Control Branch of the Office of Advanced Research and Technology, 
was used in programming the linear equations for the URIVAC 1108 digital 
computer. For the axis parallel to b^, the initial attitude angle rate 
9^ is 10 radian per second. This rate is considerably larger than those 
expected to be encountered from disturbance torques for a deep space 
mission. All other variables, i.e., 0., , etc., have an initial value 

of zero. Figures 28-30 show the attitude behavior for cases when 
i = 1, 2, 3 is in the stable range stated earlier. In particular, 

“si * “S2 - “S3 ■ 

At approximately 100 seconds, 0^ reaches its maximum amplitude of 6.4 
X 10 ^ radian. For a time of 2800 seconds, the amplitude has diminished 

-4 

to nearly 7 ^ 10 radian or one-tenth the maximum aplitude. For 9g and 
0^, the response curves are of a lower frequency than for 9^^, This is 
probably a reflection of the small inertia in comparison with Ig or 
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values for Kg^, ^33 Figures 31-33 are 50, 350^ and 350 

respectively. These fall outside of the stable ranges. The 0^ responses 
are then expected to be unstable. This is seen to be the case. However, 
the curves are not reliable as true indicators of the behavior of 0^ since 
the linear equations are not really valid for large 0^. 0^ is indeed 

unstable, but its time history probably departs markedly from the curves 
of Figures 31-33. 


Conclusions ; 

The work presented here has demonstrated the feasibility of using electric 
propulsion engines for attitude control as well as for propulsion during 
the low-thrust cruise phase of a deep space mission. In addition, it has 
been shown that flexibility associated with large roll-out solar arrays 
poses no great threat to the success of such missions. Furthermore, the 
detrimental effects from torques due to "tail-wags-dog" terms can be offset. 

Structural frequencies of the order of 0.025 cycles per second were 
encountered for the solar arrays. Only for the axis of maximum inertia, 
the one parallel to b^, did the structural influence have a pronounced 


influence on the rigid body stability characteristics. There, the permis- 
sable operating range for Kg^ reduced from I6 < < 276 to 15 < Kg2 < 201, 

Little or no change in the rigid spacecraft stability ranges for Kgj^ and 
Kg^ was encountered from structural effects. This was the case despite 
the fact that for the b^ axis, two modes of flexible motion were activated. 

For the unstable cases, the 0^ response seemed to have a growth rate 
far in excess of the 0^ or 0^ responses. A possible explanation for this 
is that Ig is greater than either or and consequently, the control 
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system is attempting to handle a large inertia torque Actually, 

is very nearly the same as However, in addition to the inertia 

torque, the control system must contend with the “tail -wags -dog" effect 
which is absent in the dynamics. The b^ axis unstable response seems 

to be intermediate to the b^ and b^ axes responses. Even though is quite 
small when compared to or there is a "tail-wags-dog" effect for this 
axis. Thus, it appears that outside of the stable operating ranges, the 
"tail-wags 'dog" influence is felt very strongly. 
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